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ABSTRACT 

Integral field spectroscopy of galaxies at redshift z ~ 2 has revealed a population of early-forming, 
rotationally-supported disks. These high-redshift systems provide a potentially important clue to the 
formation processes that build disk galaxies in the universe. A particul arly well-studied example is the 
z = 2.38 galaxy BzK— 15504, which was shown bv lGenzel et al.l (|2QQ6l ) to be a rotationally supported 
disk despite the fact that its high star formation rate and short gas consumption timescale require a 
very rapid acquisition of mass. Previous kinematical analyses have suggested that z ~ 2 disk galaxies 
like BzK— 15504 did not form through mergers because their line-of-sight velocity fields display low 
levels of asymmetry. We perform the same kinematical analysis on a set of simulated disk galaxies 
formed in gas-rich mergers of the type that may be common at high redshift, and show that the 
remnant disks display low velocity field asymmetry and satisfy the criteria that have been used to 
classify high-redshift galaxies as disks observationally. Further, we compare one of our remnants to 
the bulk properties of BzK— 15504 and show that it has a star formation rate, gas surface density, and 
a circular velocity-to- velocity dispersion ratio that matches BzK— 15504 remarkably well. We suggest 
that observations of high-redshift disk galaxies like BzK— 15504 are consistent with the hypothesis 
that gas-rich mergers play an important role in disk formation at high redshift. 

Subject headings: galaxies: formation - galaxies:high redshift - galaxies:kinematics and dynamics 



1. INTRODUCTION 

Understanding the formation of disk galaxies re- 
mains a primary but elusive goal in cosmology. The 
standard scenario for disk formation entails the qui- 
escent, dissipational collapse of rotating gas clouds 
within virialized dark m a tter halos dWhite fc Reesl Il978l : 
iFall fc Efstathioul Il980t iBlumenthal et all Il984ft . and 
this picture successfully explains the angular momen- 
tum content, size, and kinematica l structure of ob- 
served systems (e.g.. iMo et~aTlll998f ). However, CDM- 
based cosmological simulations of disk galaxy formation 
have had difficulty producing h igh angular momentum 
disks without large bulges (e.g., I Navarro fc White! fl993 : 
iNavarro fc Steinm etz 2000), almost certainly in part due 
to the fact that mergers a re common in Cold D ark Mat- 
ter (CDM) cosmologies dStewart et al.l I2QQ8L and ref- 
erences therein). More recent simulations have faired 
better by including strong feedback, which suppresses 
star formation and thereby reduces the collis ionless heat- 
ing a ssociated with stellar mer gers ("e.g.. | Abadi et all 
|2003USommer-Larsen et~^l2003nR,obertson et al.1 2004; 
iGovernato et alJl2QQ4( l. 

Motivated by these results, iRobertson et al.l ([2006a, 
hereafter, R06a) presented a supplemental "merger- 
driven" scenario for the cosmological disk galaxy forma- 
tion where extremely gas-rich mergers at high-redshift 
lead to rapidly-rotating remnant systems wit h large 



gaseou s and stellar disks. Buildi ng on work bv iBarnes 
(j2002h . iBrook et al.1 (|2QQ4f ). and iSpringel fc Hernquist 
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(|2QQ5h . R06a used a suite of hydro dynamical simulations 
to study gas rich mergers of various gas fractions, mass 
ratios, and orbital properties, and concluded that hier- 
archical structure formation and the regulation of star 
formation may work in concert to build disk galaxies 
through gas-rich mergers. Subsequent cosmological sim- 
ulations have reported the format ion of Milky Way-lik e 
disk galaxies in gas- r ich m ergers ([Governato et al.l l2007). 
and iHopkins et al.l ((2008) have examined merger rem- 
nants in order to gain a phenomenological understand- 
ing of disk survival in gas-rich mergers. Given that high 
gas fractions are the essential ingredient in this scenario, 
observations at high-redshift (when gas fractions were 
higher and mergers more common) provide an important 
testing ground for models of disk galaxy formation. 

Only recently have detai led observation s of high- 
redshift disks been possible. lErb et aT] (|2003l ) measured 
the rotation curves of UV-selected z ~ 2 galaxies (the 
BM/ BX sample, lAdelberger et al.l 120041 : ISteidel et ail 
|2004f) with YLa slit spectroscopy ( see also lErb et all 
I2006j ). iForster Schreiber et al.1 (120061) used the SINFONI 



integral field spectrograph (|Eisenhauer et al.l l2QQ3h to 
measure spatially-resolved kinematics of 14 BM/BX 
galaxies, and found rotating systems with large spe- 
cific angular momenta but significant random motions 
(u/<7~2-4). 

Among the best-studied high-redhift disks is 
BzK— 15504 at z = 2.38, which was found using 
the B z K photometric select ion technique (|Daddi et al.l 
2004). iGenzel et all (|2006l . hereafter G06) used the 
SINFONI spectrograph with adaptive optics on the Very 
Large Telescope (VLT) to measure the kinematic struc- 
ture BzK— 15504 with an incredible spatial resolution 
of 0.15" (~ 1.2 kpc at 2; ~ 2.4). These observations re- 



vealed massive (M* - 8 x 10 10 M o , M gas - 4 x 10 10 M o ) 
rotating disk galaxy with a large star formation rate 
(SFR ~ 14OM yr -1 ) and substantial random motions 
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(v/cr ~ 3). The inferred gas consumption timescale for 
BzK— 15504 is only ~ 285 Myr, or about one-tenth the 
age of the universe at z ~ 2.38, and implies a remarkably 
rapid acquisition of its mass. Nonetheless G06 conclude 
that BzK— 15504 is not a merger remnant because it 
displays velocity fields similar to those expected for 
quiescent disks. 

The existence of high-redshift disks that have not expe- 
rienced large mergers may be surprising in light of ACDM 
expectations. The merger rate per dark matter halo is 
predicted to be a factor of ~ 15 times higher at z ~ 2.4 
than at z ~ (jFakhouri fc Mai l2Q08h . Approximately 
10% of galaxy-size halos at z ~ 2.4 should have under- 
gone a nearly one-to-one merger ( > 1 : 1.25) in the last 
~ 500 Myr, and ~ 50% should have experienced a > 1 : 3 
merger in the last ~ 1 Gyr (Stewart, Bullock et al., in 
preparation). However, these mergers are expected to be 
gas rich. Observations of galaxies at z ~ 2 imply gas 
fractions ~ 10 — 20 times higher than for nearby galax- 
ies at fixed stellar mass (|Erb et al.l I2QQ6I : iGavazzi et al.l 
2008). These two pieces of information motivate us to 
consider gas-rich mergers as a means to explain z ~ 2 
galaxies with young stellar populations, short gas ex- 
haustion timescales, and disk-like kinematics. 

Any model that attempts to explain the properties of 
high-redshift disks must not only reproduce their bulk 
properties (stellar masses, rotational velocities, gas frac- 
tions, star formation rates, etc.) but also their detailed 
kinematic properties. An important kinematic metric for 
cl assifying high-redshi ft galaxies as disks was developed 
bv lShapiro et al.l (|2008l he reafter S08), who exten ded the 
"kinemetry" technique of iKrajnovic et ahl (|2006l ) in or- 
der to examine the velocity fi elds of z ~ 2 galaxies from 
IForster Schreiber et all (j2QQ6h and G06. Specifically, S08 
provided a straightforward means to classify z ~ 2 galax- 
ies as disks based on the symmetry of their velocity fields. 

The purpose of this Letter is to determine whether 
galaxies formed in a gas-rich mergers can match the ob- 
served properties of z ~ 2 disks. We describe our simu- 
lations in ^2] and the kinematical analysis in In ^Hwe 
demonstrate that a set of gas-rich merger remnants sat- 
isfy the disk galaxy kinematical classification developed 
by S08, and compare one of the remnants in particular to 
the detailed properties of BzK— 15504. Overall, we find 
good agreement between the properties of simulated disk 
galaxies formed in gas-rich mergers and observed z ~ 2 
disks like BzK— 15504. We assume a ACDM cosmology 
with Q m ~ 0.3, ~ 0.7, and a Hubble parameter of 
H = TOkms^Mpc" 1 . 

2. SIMULATION METHODOLOGY 

The simulations studied here are culled from the 
simulation suite presented in R06a. For details 
of the simulation meth odolog y we refer the reader 
to ISpringel fe Hernquistl (l2QQ3h. ISpringel et all (|2005f ). 
ISpringell (|2005f h and iRobertson et al.l (|2006bl jd). but a 
summary follows. The simulations are calculated using 
the N-body and Smoothed Particl e Hydrodynam ics code 
GADGET2 (jSpringel et aHl200ll : ISpringelll2005h . Each 
binary merger occurs during a parabolic encounter be- 
tween otherwise isolated disk galaxies. The progenitor 
galaxy models contain gaseous and stellar disks embed- 
ded in massive dark matter halos, with structural pa- 
rameters scaled appropriately for disk galaxies in the 



ACD M cosmology (e.g., iMo et al.l fl998b ISpringel etafl 
l2QQ5h . The gas and stars are initialized in exponential 
disks with a scale height Hd to scale length Rd ratio 
of Hd/Rd — 0.2, while the dark matte r halo is initial- 
ized as rotating (A = 0.033) Hernquist (1990) distribu- 
tion with an approximate iNavarro et al] (|l996l ) concen- 
tration of cnfw = ^200 / Rs = 9 (the simulations of R06a 
were initialized for concentrations typical of z ~ 0). We 
note that while the concentration is higher tha n that ex- 
pecte d for similar mass halos at z ~ 2 (e.g., Bull ock et al.l 
l200ll ). the integrated mass within the central ~ 10 kpc 
is approximately correct. The evolution in halo con- 
centrations with redshift is primarily driven by growth 
in halo virial radii with time, while halo central densi- 
ti es remain approxima tely constant (see, e.g., Figure 18 
of IWechsler et al.ll2002j ). The progenitor galaxy models 
each contain 40,000 gas, 40,000 stellar disk, and 180,000 
dark matter halo particles, with a gravitational softening 

of e = lOOft-V- 

The prescription for star formation and interstellar 
medium (ISM) physi c s in the simulations follows the 
ISpringel fc Hernquistl (|2003l ) model that implements a 
sub-resolution treatment of supernova feedback and the 
multiphase ISM. Star formation operates in the dense 
ISM, following a ISchmidtl (|l959f ) law where the star 
formation timescale decreases with the local dynamical 
time. The simulations of R06a utilize the strong su- 
perno va feedback supplied by the ISpringel fc Hernquistl 
(|2003l ) multiphase ISM model to suppress star forma- 
tion sufficiently during a merger to allow remnant disks 
to form. The simulations do not include an ultraviolet 
(UV) background, but we note that the increased UV 
background at z ~ 2 would tend to suppress the overall 
gas consumption during the merger and may act to in- 
crease the rotational support of the final remnants. For 
comparison with the observations of BzK— 15504 by G06, 
we focus on a single equal-mass coplanar binary merger of 
disk galaxies with virial velocities V v { r = 160kms _1 and 
a pericentric passage distance of r per i ~ 2Rd = 8.6kpc. 
The model is simulation " GE" in the study by R06a (see 
their Table 1) and is s i milar to the simulation studied by 
ISpringel fc Hernquistl (|2005f ) . At the time of final coales- 
cence At « 450 Myr after the first passage, the merging 
system is « 60% gas and experiences a large burst of 
star formation. The system averages SFR w 15OM0yr -1 
until reaching a quiescent rate of SFR « 15M yr _1 
approximately 285 Myr later. We analyze the simula- 
tion on the declining peak of the vigorous star formation 
burst, approximately 140 Myr after the start of the fi- 
nal coalescence (roughly 100 Myr before the third panel 
of Figure 1 in R06a). By this stage, a large gaseous 
disk has formed from the residual angular momentum of 
the merger. Eventually, the system forms a rotationally- 
supported stellar component with a majority of its stellar 
mass distributed in thin and thick disks (see Figures 2 
and 3 and Table 2 of R06a). 

3. KINEMATICAL ANALYSIS 

The kinematical properties of the simulated merger 
remnants are analyzed using a method designed to 
closel y approximate the integ r al fiel d spectroscopy stud- 
ies of IForster Schreiber et all (|2006l ) and G06. Figure [1] 
illustrates this analysis for the GE remnant from R06a 
discussed in §2. The disk is inclined to i = 48deg and 
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Fig. 1. — Density and velocity fields of a simulated disk galaxy formed in a gas-rich merger 
i = 48 deg and a position angle PA = 24 deg West of North for comparison with the redshift z 



The disk is rotated to an inclination of 
2.348 galaxy BzK-15504 (G06, Figure 

3). The scale bar in the upper panels shows the relative size of 0.5" at z = 2.348. The left, middle, and right columns show the surface 
mass density (SMD), velocity (v\ os ), and velocity dispersion (<r los ), respectively, for gas associated Ha: emission. The color scales display 
the range of SMD, velocity, and velocity dispersion for all rows. The upper row (a-c) provides high-resolution maps of the simulation, 
while the middle row (d-f) is smoothed using a 0.15" FWHM gaussian, typical of VLT resolution with adaptive optics. The bottom row 
(g-i) shows the sm oothed fields binned to a 0.05" pixel scale, chosen to match the capabilities of the SINFONI integral field spectrograph 
(Eisenhauer et al. 2003). Gaussian noi se is a dded to the binned maps to match the typical velocity errors reported by S08. The kinemetric 
analysis designed by Krainovic et al. (2006) and S08 is performed on the pixelated i;i os (panel h) and cri os (panel i) maps and yields an 
asymmetry parameter of K asym = 0.1. S08 classify galaxies with K asym < 0.5 as disks. 

rotated to a major kinematic axis position angle of PA = 
24 deg West of North in order to match BzK— 15504 (see 
G06, their Figure 3). Star- forming gas and diffuse gas 
with a temperature in the range 3 x 10 3 < T < 3 x 10 4 K 
are selected to approximate the observed Ha-emitting 
gas and the size of the simulated galaxy is scaled to the 
redshift of BzK— 15504, z = 2.38 (see the angular scale 
bar in panels a-c of Figure [T]). The three columns of Fig- 
ured] show gas surface mass density (SMD), line-of-sight 
velocity v\ os , and velocity dispersion <j\ os maps of the 
simulated disk gas displayed at three different smooth- 
ing scales. The upper row (a, b, c) is rendered at high 
resolution. Note the clear disk morphology (a) and disk- 
like velocity field (b), while the random motions (c) are 
larger than for thin disks in the local universe (we mea- 
sure, v\ os /<j\ os ~ 2 — 4 in the remnant gas disk). In the 
middle row (d, e, f) the fields are smoothed by a gaussian 
kernel with a FWHM= 0.15" in order mimic the angular 
resolution obtained by G06. 5 In the lower row (g, h, i) 
the smoothed maps are then binned to the 0.05" pixel 
scale of the SINFONI observations by G06 and gaussian 



noise is added to match the typical velocity errors re- 
ported by S08 (see their Figure 6). 

The kinemetry analysis developed bv iKrajnovic et all 
(2006)) is then applied to the pixelated ^i os and ai os maps. 
This method calculates an expansion of the velocity field 
along an ellipse with semi-major axis a as a function of 
azimuthal angle i/j of the form 



K(a^) = Aq (a) + [A n (a) sin nijj + B n (a) cos nip] , 



5 Note that most kinematical studi es of z ~ 2 galaxies have 
' 0.5" resolution (see Forster Schreiber et al. 2006 and S08). 



(i) 

where n max is a maximum term in the expansion and 
the coefficient B\ describes the circular velocity of the 
system. For high-quality kinematic al data of nearby 
galaxies (e.g., from SAURON, see iBacon et all 120011 : 
ICappellari et al.l l2006). the kinemetry analysis provides a 
detai led exposition of the k inematical structure of galax- 
ies ([Krajnovic et al.l 12006). The comparatively coarser 
data available for z ~ 2 galaxies motivated S08 to use 
a mean kinemetric analysis where the asymmetries are 
characterized by averages of the amplitude coefficients 
in Equation [TJ Re- writing the kinemetry coefficients as 
k n = y / A^ + S08 introduced average velocity and 
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Properties of BzK 



TABLE 1 

-15504 vs. Simulated Disk Merger Remnant 



Property 1 


BzK-15504 1 


Disk Remnant 2 


Dynamical Mass (r< 1.1")' 


1.1 ± 1 x 10 n M Q 


1.2 x 10 n M Q 


Stellar Mass 


7.7(+3.9, -1.3) x 1O 1O M 


7.3 x 1O 1O M 


Gas Mass 


4.3 x 1O 1O M 


3.8 x 1O 1O M 


Average E gas 


35OM pc- 2 


3O6M pc- 2 


Average £ S fr 


l^MQyr-^pc -2 


l.lMQyr-^pc -2 


Vc/cr 


3±1 


3.2 ±2 


Gas Disk Scale Length 


4.5±lkpc 


2.6 kpc 



References. — (1) IGenzel etaLl (120061 ); (2) IRobertson et al.l (|2006a) 
velocity dispersion asymmetry parameters 

__ / ^n=2 ^n,v/4 \ _ / E n =l ^n,a/5 \ 



(2) 

where the subscripts v and a refer to quantities calcu- 
lated from the v\ os and cri os maps, and the subscript r in- 
dicates that the average is performed over the ellipses at 
different radii used in the kinemetry analysis. We use the 
kinemetry code made available by D. Krajnovic to calcu- 
late v asym and cr asym for the simulated disk remnant. The 
kinemetry is performed over six ellipses separated by ~ 3 
pixels in the map, roughly corresponding to the resolu- 
tion of the s imulated image to minimiz e beam smearing 
effects (e.g.. Ivan den Bosch et al.ll2Q00h . The ellipse po- 
sition angles and aspect ratios are fixed to reflect the best 
fit values for BzK— 15504 adopted from G06. The num- 
ber, spacing, and orientation of the ellipses can change 
the precise values of v asym and cr asym , but do not alter 
the conclusions or interpretations of this work. Increas- 
ing the level of noise in the simulated maps will tend to 
increase the values of v asym and cr asym , and would likely 
degrade the capability of the simulations to match the 
observations rather than improv e our results. We note 
that while the kinemetry method Kraj novic etall (2006) 
has pr eviously been used analyze simulated merge r rem- 
nants (Jessei tet al.ll2007t iKronberger et aT]|2007t ). only 
very gas-rich mergers are expected to be associated with 
disk formation and kinemetry has not yet been performed 
on such systems. 

4. RESULTS AND DISCUSSION 

The bottom row of Figure [1] shows the SMD, i?i os , and 
<j\ os maps used to perform the kinematical analysis of the 
remnant disk. The average kinemetric parameters v asym 
and cr asym (Equation 2) are measured from the velocity 
maps. For the disk merger remnant plotted in Figure 
Q] we find v asym = 0.076 and cr asym = 0.063, which re- 
flect the relative symmetry the velocity fields. The com- 



bined asymmetry parameter K asym , asym , asym 

provides a global characterization of asymmetry in the 
kinematical structure of the galaxy, and S08 suggest 
K asym < 0.5 as an observational criterion for distinguish- 
ing between disk galaxies and mergers. The disk rem- 
nant examined in Figure [T] satisfies this criterion with 
K asym = 0.1, even though it was formed in an equal 
mass, gas-rich merger. We have performed the same 
analysis on a range of other simulated merger remnant 
disk galaxies from R06a. In each case, we examine the 



merger product soon after the merger and 150 Myr after 
the time of the final coalescence. We find similar results 
for disk remnants formed in gas-rich major polar mergers 
(R06a simulation GCoF, K asym = 0.14), major inclined 
mergers (R06a simulation GCoO, K asym = 0.15), and 
minor (1:8 mass ratio) coplanar mergers (R06a simula- 
tion FCm, if asym = 0.16). If real galaxy analogues of 
these simulated galaxies were observed with SINFONI, 
a kinematical analysis would likely classify them (cor- 
rectly) as disks. We emphasize that in each case the gas 
disk develops rapidly, with v/cr and K asym both declining 
by factors of ~ 2 — 3 within ~ 100 Myr of the height of 
merger coalescence. 

While both the simulated disk merger remnants and 
observed z ~ 2 disk galaxies satisfy the observational cri- 
terion of S08 designed to identify disk galaxies, the sim- 
ulations should also closely match other observed prop- 
erties of z ~ 2 disks. Consider the galaxy BzK— 15504 
observed by G06, which has an inferred star formation 
rate of SFR « 14OlgJ°M yr _1 and a dynamical mass of 
M dyn = 1.1 zb 0.1 x 10 n M o within a radius of r < 9kpc. 
The simulated disk remnant analyzed in Figure [T] has 
an average star formation rate of SFR w 15OM yr =1 
over the previous 100 Myr and a dynamical mass of 
M dyn = 1.2 x 10 n M o . The inferred stellar and gas 
masses for BzK-15504 are M* = 7.7±?;jj x 10 10 M o and 
M gas = 4.3 x 1O 1O M , while the disk remnant has a 
stellar mass of M* = 7.3 x 1O 1O M and a gas mass of 
M gas = 3.8 x 1O 1O M . Table CD compares the salient 
properties of BzK— 15504 and the simulated disk rem- 
nant, and they are similar in each case. Overall, the 
galaxy BzK— 15504 displays properties that are remark- 
ably similar to the example simulated disk galaxy merger 
remnant shown in Figure [TJ 

We note that the properties of z ~ 2 disks like 
BzK— 15504 are rather extreme compared with local, iso- 
lated disk galaxies. The gas surface density, star for- 
mation rate, and random motions (v/cr) of BzK— 15504 
are roughly an order of magnitude larger than for the 
Milky Way, even as their rotational velocities are sim- 
ilar a nd their stellar masses differ by less than 50% 
(e.g.. iBlitz et al.1 I l980l: iRanal Il99ll: iKent et all Il99ll : 



iKlvpin et all 120021 : iLevine et al. 2006). Our simulated 
disk remnant displays similar properties (see Table [1]) 
owing to the dynamical effects of the merger from which 
it formed. It is possible that such a pre-formed (hot) disk 
would tend to settle into a more classical Milky- Way type 
thin disk at late times if the pre-formed stars adjust to 
to the infall of (quiescent) gas-disk material, as might 
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be expected in a more classical stage of disk formation. 
Subsequent gaseous infall from the intergalactic medium 
may also help prolong star formation in the remnant and 
increase the gas disk scale length. 

While the observed properties of z ~ 2 disks resem- 
ble simulated disk systems formed in gas-rich mergers, 
other viable exp l anations for their origin exist. Notably, 
iBournaud et all (|2007l |2008[ ) have suggested that high- 
redshift disk galaxies undergo a "clump-cluster" phase 
where large-scale gravitational instability leads to a tur- 
bulent, cl umped gaseous distribution (see also iNoguchil 
Il998l . fl999l ). An analysis of clump- cluster ga l axy m odels 
using kinemetry method of iKrajnovic et all ([2006) and 
S08 would provide a useful test of this scenario. 

As the quantity and quality of kinematical data for 
z ~ 2 disk galaxies improve, the observational con- 



straints will better distinguish between various models 
for their formation. These observations will be essential 
for determining how disk galaxies are assembled over the 
history of cosmic structure formation. 



BER gratefully acknowledges support from a Spitzer 
Fellowship through a NASA grant administrated by the 
Spitzer Science Center, and the Kavli Institute for Cos- 
mological Physics at the University of Chicago. JSB was 
supported by NSF grants and the Center for Cosmol- 
ogy at UC Irvine. We also thank D. Krajnovic for mak- 
ing his kinemetry analysis code easily accessible and A. 
V. Kravtsov for useful comments and suggestions. We 
appreciate helpful comments and suggestions from the 
anonymous referee. 



REFERENCES 



Abadi, M. G., Navarro, J. F., Steinmetz, M., & Eke, V. R. 2003, 
ApJ, 591, 499 

Adelberger, K. L., Steidel, C. C, Shapley, A. E., Hunt, M. P., Erb, 

D. K., Reddy, N. A., & Pettini, M. 2004, ApJ, 607, 226 
Bacon, R., et al. 2001, MNRAS, 326, 23 
Barnes, J. E. 2002, MNRAS, 333, 481 

Blitz, L., Fich, M., & Stark, A. A. 1980, in IAU Symposium, Vol. 87, 
Interstellar Molecules, ed. B. H. Andrew, 213-217 

Blumenthal, G. R., Faber, S. M., Primack, J. R., & Rees, M. J. 
1984, Nature, 311, 517 

Bournaud, F., et al. 2008, ArXiv e-prints, 0803.3831 

Bournaud, F., Elmegreen, B. G., & Elmegreen, D. M. 2007, ApJ, 
670, 237 

Brook, C. B., Kawata, D., Gibson, B. K., & Freeman, K. C. 2004, 

ApJ, 612, 894 
Bullock, J. S., et al. 2001, MNRAS, 321, 559 
Cappellari, M., et al. 2006, MNRAS, 366, 1126 
Daddi, E., Cimatti, A., Renzini, A., Fontana, A., Mignoli, M., 

Pozzetti, L., Tozzi, P., & Zamorani, G. 2004, ApJ, 617, 746 
Eisenhauer, F., et al. 2003, in Proceedings of the SPIE, Vol. 4841, 

1548-1561 

Erb, D. K., Shapley, A. E., Steidel, C. C, Pettini, M., Adelberger, 
K. L., Hunt, M. P., Moorwood, A. F. M., & Cuby, J.-G. 2003, 
ApJ, 591, 101 

Erb, D. K., Steidel, C. C, Shapley, A. E., Pettini, M., Reddy, 

N. A., & Adelberger, K. L. 2006, ApJ, 646, 107 
Fakhouri, O., & Ma, C.-P. 2008, MNRAS, 386, 577 
Fall, S. M., & Efstathiou, G. 1980, MNRAS, 193, 189 
Forster Schreiber, N. M., et al. 2006, ApJ, 645, 1062 
Gavazzi, G., et al. 2008, A&A, 482, 43 
Genzel, R., et al. 2006, Nature, 442, 786 
Governato, F., et al. 2004, ApJ, 607, 688 

Governato, F., Willman, B., Mayer, L., Brooks, A., Stinson, G., 
Valenzuela, O., Wadsley, J., & Quinn, T. 2007, MNRAS, 374, 
1479 

Hernquist, L. 1990, ApJ, 356, 359 

Hopkins, P. F., Cox, T. J., Younger, J. D., & Hernquist, L. 2008, 

ArXiv e-prints, 806 
Jesseit, R., Naab, T., Peletier, R. F., & Burkert, A. 2007, MNRAS, 

376, 997 

Kent, S. M., Dame, T. M., & Fazio, G. 1991, ApJ, 378, 131 
Klypin, A., Zhao, H., & Somerville, R. S. 2002, ApJ, 573, 597 



Krajnovic, D., Cappellari, M., de Zeeuw, P. T., & Copin, Y. 2006, 

MNRAS, 366, 787 
Kronberger, T., Kapferer, W., Schindler, S., & Ziegler, B. L. 2007, 

A&A, 473, 761 
Levine, E. S., Blitz, L., & Heiles, C. 2006, ApJ, 643, 881 
Mo, H. J., Mao, S., & White, S. D. M. 1998, MNRAS, 295, 319 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 

563 

Navarro, J. F., & Steinmetz, M. 2000, ApJ, 538, 477 
Navarro, J. F., & White, S. D. M. 1994, MNRAS, 267, 401 
Noguchi, M. 1998, Nature, 392, 253 
Noguchi, M. 1999, ApJ, 514, 77 
Rana, N. C. 1991, ARA&A, 29, 129 

Robertson, B., Bullock, J. S., Cox, T. J., Di Matteo, T., Hernquist, 

L., Springel, V., & Yoshida, N. 2006a, ApJ, 645, 986 
Robertson, B., Cox, T. J., Hernquist, L., Franx, M., Hopkins, P. F., 

Martini, P., & Springel, V. 2006b, ApJ, 641, 21 
Robertson, B., Hernquist, L., Cox, T. J., Di Matteo, T., Hopkins, 

P. F., Martini, P., & Springel, V. 2006c, ApJ, 641, 90 
Robertson, B., Yoshida, N., Springel, V., & Hernquist, L. 2004, 

ApJ, 606, 32 
Schmidt, M. 1959, ApJ, 129, 243 

Shapiro, K. L., et al. 2008, ArXiv e-prints, 0802.0879 
Sommer-Larsen, J., Gotz, M., & Portinari, L. 2003, ApJ, 596, 47 
Springel, V. 2005, MNRAS, 364, 1105 

Springel, V., Di Matteo, T., & Hernquist, L. 2005, MNRAS, 361, 
776 

Springel, V., & Hernquist, L. 2003, MNRAS, 339, 289 
— . 2005, ApJ, 622, L9 

Springel, V., Yoshida, N., & White, S. D. M. 2001, New Astronomy, 
6, 79 

Steidel, C. C, Shapley, A. E., Pettini, M., Adelberger, K. L., Erb, 

D. K., Reddy, N. A., & Hunt, M. P. 2004, ApJ, 604, 534 
Stewart, K. R., Bullock, J. S., Wechsler, R. H., Mailer, A. H., & 

Zentner, A. R. 2008, ApJ, Accepted, 711 
van den Bosch, F. C, Robertson, B. E., Dalcanton, J. J., Sz de 

Blok, W. J. G. 2000, AJ, 119, 1579 
Wechsler, R. H., Bullock, J. S., Primack, J. R., Kravtsov, A. V. & 

Dekel, A. 2002, ApJ, 568, 52 
White, S. D. M., & Rees, M. J. 1978, MNRAS, 183, 341 



